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Abstract. The Large Area Telescope (LAT) on 
the Fermi Gamma-ray Space Telescope (formerly 
GLAST), with its improved sensitivity relative to 
previous generation gamma-ray telescopes, is signif- 
icantly increasing the number of known gamma-ray 
sources in the sky, including pulsars. In addition to 
searching for gamma-ray emission from known radio 
pulsars, it is now possible to successfully perform 
blind searches for pulsars on the gamma-ray data 
alone, with the goal of uncovering a new population 
of potentially radio-quiet (Geminga-Uke) pulsars. We 
describe our methods and some recent results from 
our searches, performed using the time-differencing 
technique, which have resulted in the discovery of a 
large number of new gamma-ray pulsars. 
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I. Introduction 

Prior to the launch of Fermi, only six pulsars 
were known to emit > 100 MeV gamma rays. The 
high energy pulsations from these spinning neutron 
stars were all detected by folding the gamma-ray data 
with the known ephemeris from radio observations, 
or in the case of the radio-quiet Geminga, using the 
X-ray ephemeris. In addition to these six gamma-ray 
pulsars, the EGRET mission on the Compton Gamma 
Ray Observatory (1991-2000) detected hundreds of 
sources, a majority of which remained unidentified |1|. 
Many unidentified gamma-ray sources, in particular 
those in the Galactic plane, have long been suspected 
of being pulsars, but deep radio and X-ray searches 
for pulsations have often come up empty. One of 
the key scientific objectives of Fermi is to resolve 
the gamma-ray sky and identify these sources. In the 
first months of the mission, Fermi is meeting those 
expectations and uncovering an entire population of 
pulsars hitherto unknown at any wavelength. 

II. The Fermi Large Area Telescope (LAT) 

The Fermi satellite, consisting of the Gamma-ray 
Burst Monitor (GBM) and the Large Area Telescope 
(LAT), was launched on 11 June 2008 into a low Earth 
circular orbit at an altitude of 565 km and an inclination 
of 25.6°. The LAT Q is a pair-production telescope 
with large effective area (^8000 cm^) and field of view 
(2.4 sr), sensitive to gamma rays between 20 MeV 
and > 300 GeV. Although its commissioning phase 
(30 June to 30 July 2008) was primarily intended for 



instrument checkout and calibration, several scientific 
results were obtained with these early data, including 
the first discovery of a radio-quiet gamma-ray pulsar 
in a blind search |3|. The LAT began nominal science 
operations on 11 August 2008, and has since been 
observing mostly in survey mode, scanning the entire 
sky every three hours. The overall sensitivity of the 
LAT is '^25 times that of EGRET, while the angular 
resolution is also significantly improved (it ranges from 
-3-6° at 100 MeV to -0.1-0.2° at 10 GeV). The 
mission was designed with a five-year lifetime and a goal 
of at least ten years of operations. The scientific goals of 
the mission span a wide range of topics, including un- 
derstanding the emission mechanisms in Active Galactic 
Nuclei (e.g. |4|, |5|, |6|), pulsars (e.g. Q, HI, |l9j), and 
supernova remnants (SNR) (e.g. QO), fH, (Til), the 
high energy emission from Gamma-ray Bursts (e.g. lfT3l . 
El, CSl), and the nature of dai'k matter (e.g. HH, ifTTl . 
CD. 

III. Blind Searches for Gamma-Ray Pulsars 

Detecting gamma-ray pulsars is challenging for a 
variety of reasons. The low count rates require long 
integration times (up to several years) and result in very 
sparse data sets. This, combined with the short spin 
periods of most pulsars (<ls), presents serious com- 
putational problems for standard fully coherent search 
techniques involving Fast Fourier Transforms (FFT). In 
addition, during such long time spans, pulsars, especially 
young ones, will often experience irregularities in their 
timing behavior, such as "timing noise" or glitches (e.g. 
PSR B1706-44 further complicating the problem. 

The term "blind" refers to the fact that the spin 
parameters of the potential pulsar (e.g. frequency, /, 
and frequency derivative, /) are unknowrQ. A relatively 
small number of locations in the sky are searched 
for pulsations. These include promising sources from 
multiwavelength observations, as well as the locations 
of LAT-detected gamma-ray sources. A list of the most 
significant LAT sources detected in the first three months 
(> lOcr) was recently released to the public ll20l . 
All those not associated with AGN were searched. In 
addition, we searched other less significant LAT sources, 
from an internal catalog, which we deemed interesting 
pulsar candidates. At least one pulsar was found from a 
LAT source not included in the bright source list II2TII . 

The arrival times of photons must be translated to the 
Solar System barycenter assuming a location, and for 

' Unlike in radio searches, we need not be concerned with searching 
over a range of values of dispersion measure (DM). 
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Fig. 1. LAT gamma-ray image, in celestial coordinates, of the Cygnus 
region of the Galactic plane, showing three bright sources, coincident 
with formerly unidentified EGRET sources 3EG J2021+3716, 3EG 
J2020+4017, and 3EG J2033+4118. All three are new gamma-ray 
pulsars : the known radio pulsar PSR J2()2 1 -1-365 1 1 22 1 and two gamma- 
ray pulsars discovered in the blind search, LAT PSR J2()2 1+4044 in 
the Gamma Cygni SNR and LAT PSR J2032-H4127 in the Cygnus 
OB2 association 1211 . 



each source, a large parameter space in / and / must 
be covered. Using the Crab to define the upper limit 
for young pulsars, we chose the maximum frequency of 
64 Hz and a minimum ratio of ///=-L25 x 10^^^ s^^ 
(the parameter space is explored in steps of ///). These 
parameters cover roughly 85% of the ~2000 pulsars 
contained in the ATNF database. 

Figure 1 shows a LAT gamma-ray image of the 
Cygnus region of the Galactic plane, and serves to 
illustrate the different ways in which the LAT is finding 
gamma-ray pulsars. The bright source at the bottom is 
the gamma-ray counterpart of the known radio pulsar 
PSR J2021H-3651 (also known as the "Dragonfly"). This 
was the first gamma-ray pulsar detected by the LAT f22l, 
and indeed the first new "post-EGRET" gamma-ray 
pulsar, although the original discovery was made slightly 
earlier by AGILE ll23l . As in the case of the six 
EGRET pulsars, the gamma-ray pulsations were detected 
by folding the gamma-ray photons using the known 
radio ephemeris of the pulsar. Above the "Dragonfly" 
pulsar lies LAT PSR J2021H-4044 in the Gamma Cygni 
SNR. This source had long been suspected of being 
a pulsar and deep Chandra observations led to a list 
of possible X-ray counterparts [24 1, which we used to 
search for, and find, the gamma-ray pulsar. Finally, in 
the top right corner of Figure 1, LAT PSR J2032H-4127 
is a pulsar coincident with the Cygnus OB2 association. 
The pulsations in this case were discovered in a blind 
search using the location of the LAT source. These three 
pulsars, like the majority of the gamma-ray pulsars dis- 
covered in our blind searches so far, are coincident with 
formerly unidentified EGRET gamma-ray sources: 3EG 
J2021H-3716, 3EG J2020H-4017, and 3EG J2033H-4118 
respectively. Furthermore, LAT PSR J2032-I-4127 is also 
coincident with the first unidentified TeV source, discov- 
ered originally by the HEGRA experiment ll25l . 
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Fig. 2. LAT PSR J2032+4127 contours of detection significance 
over a range of period and period derivative, as obtained by running 
the "prepfold" command on PRESTO |28 |. A high value of reduced 
is shown in red and a lower value in white or purple. t2H . 
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Fig. 3. A ~5 ks Swift observation showing an X-ray image 
of the region of the sky, in celestial coordinates, around LAT PSR 
J0633+0632. The large circle represents the LAT 95% error circle for 
the source after three months, while the smaller circle is the LAT 
95% error circle after 6 months. The crosshair marks the location of 
the assumed X-ray counterpart of the pulsar, which we have named 
Swift J063344+0632.4 



A. The time-differencing technique 

As described in the previous section, the low fluxes 
typical of gamma-ray sources result in very long, sparse 
data sets. The application of traditional EFT techniques 
on such data quickly become unfeasible, for several 
reasons. First, the frequency resolution from such long 
viewing periods results in FFTs with billions (or tens 
of billions) of frequency bins. More importantly, the 
relatively large frequency derivatives of pulsars means 
that these FFTs would have to be computed a large 
number of times in order to span a realistic parameter 
space, given the small step size in / required to keep the 
signal power within a single EFT bin. To make matters 
worse, the long viewing periods increase the chances 
that the pulsar will experience a glitch during this time. 
A new technique, known as "time-differencing" ll26l . 
Il27l . has been developed, which greatly improves the 
prospects of carrying out these searches. Based on the 
premise that a periodic signal present in a time series 
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(a) Phase plot obtained using the LAT source location. (b) Phase plot obtained using the Swift source location. 

Fig. 4. LAT PSR J0633+0632 phase plots and evolution of the reduced with time, using PRESTO (281. The plots illustrate the strength 
of the signal over the course of the entire observation, and the crucial importance of a precise location in determining a good timing solution. 



should also show up when analyzing the differences of 
these times, the method involves taking time differences 
only up to a short time window (on the order of 
days, rather than months or years). In doing this, the 
required number of FFT bins, N, is greatly reduced 
(N = 2 X / X Tu], where / represents the maximum fre- 
quency searched and is the maximum time-difference 
window). The reduced frequency resolution also results 
in a larger step size required for /, greatly reducing 
the number of / trials needed. The overall result is 
that the computational and memory costs become a 
fraction of the standard FFT methods, with only a small 
reduction in sensitivity. This makes it possible to carry 
out a large number of blind searches over a realistic 
parameter space on standard desktop computers with 
only a few GB of memory. That said, the blind search 
efforts have been significantly enhanced by access to the 
UCSC Astrophysics Supercomputer Pleiades, with 207 
Dell PowerEdge 1950 compute nodes (828 processing 
cores). Despite the efficiency of the time-differencing 
technique, it is still desirable to limit the number of 
source locations in the sky that are searched, given both 
the finite computational resources available, and also to 
avoid incurring large penalties in the significance of our 
detections, due to the increasing number of trials used 
in finding a particular signal. 

After a candidate pulsar is identified in our search, 
we use standard pulsar tools, such as PRESTO and 
TEMPO to refine the ephemeris. Figure 2 shows the 
results obtained by running the prepf old command on 
PRESTO |i281, showing the significance of the detection. 
Often, the relatively large uncertainty in the LAT loca- 
tion makes it difficult to obtain a good timing solution. 



We then use several techniques to try to refine the 
location of the source, the most fruitful of which involves 
multiwavelength observations. We have been fortunate 
to obtain Swift follow-up observations of several LAT 
pulsars, leading to the identification of potential X-ray 
counterparts. Figure 3 shows the result of a short Swift 
observation of one of our newly-discovered pulsars, LAT 
PSR J0633H-0632. The largest circle represents the 95% 
error circle derived by the LAT after three months of 
data. The medium circle shows the 95% error circle after 
six months of data. Finally, the crosshair pinpoints the X- 
ray source we believe to be the counterpart to the pulsar, 
a source we have named Swift J063344-I-0632.4. Figure 
4 illustrates the crucial importance of having a precise 
location. The timing solution used in the left panel was 
the best we were able to obtain (using only a frequency 
and frequency derivative) for the LAT 6-month location. 
The panel on the right represents the best timing solution 
when using the location of the Swift X-ray source to 
barycenter the gamma-ray photons. 

IV. A NEW POPULATION OF GAMMA-RAY PULSARS 

The first major Fermi discovery was the detection of 
pulsations from CTA 1 in a blind search |3|. This young, 
nearby, shell-type SNR was discovered in radio in the 
1960s and X-ray observations show a well-localized 
central point source, RXJ0007. 0-1-7303, embedded in 
a pulsar wind nebula (PWN). High energy (> 100 
MeV) emission was detected by EGRET from 3EG 
JOO 10-1-7309, coincident with this source, but the low 
number of photons and poor localization made the search 
for pulsations very challenging. The new pulsar, with 
a period of 315.8637050 ms and a period derivative 
of 3.615x10^^'^ s s^^ is a typical young, energetic 
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Fig. 5. Folded light curve of LAT PSRJ2()32+4127, with 32 phase 
bins per period, using 5 months of data with E>300 MeV, and R<0.8° . 
Two rotations are shown for clarity f21 j. 

pulsar, with a derived characteristic age of ~ 14,000 
years (consistent with the estimated age of the SNR) 
and a spin-down power of 4.5x10'^^ erg s^^. For more 
details, see |31. Two other pulsai's, LAT PSR J 1809- 
2332 ("Taz") and LAT PSR J1836H-5925 (the "next" 
Geminga) were discovered in our blind searches dur- 
ing the commissioning phase of the instrument. All 
three of these pulsars were associated with former 
EGRET unidentified sources and in all cases it was 
long-suspected that a pulsar was responsible for the 
gamma-ray emission. Furthermore, the existence of X- 
ray sources thought to be the counterparts of such pulsars 
facilitated the detection of their gamma-ray pulsations. 
LAT PSR J1836H-5925 has been named as one of the top 
5 most interesting Galactic sources detected by Fermi 
so far 0. Our blind searches have continued to yield 
new and exciting detections of gamma-ray pulsars, to 
the point where there are now over a dozen new pulsars 
discovered through their gamma-ray pulsations |21 1. It is 
too early to say how many of these pulsars will turn out 
to be radio-quiet. Follow-up observations at the world's 
largest radio telescopes are in progress. The results of 
such searches will have important consequences in the 
population studies of pulsars and in our understanding 
of the underlying pulsar emission mechanisms. 

V. Conclusions 

The Fermi LAT has been scanning the gamma-ray 
sky with unprecedented sensitivity since its launch in 
June 2008. One of its key scientific missions is to 
resolve the gamma-ray sky by identifying gamma-ray 
sources that previous experiments left unidentified, and 
to discover new sources. The hugely improved sensitiv- 
ity and angular resolution of the LAT, combined with 
a new powerful time-differencing technique, makes it 
possible, for the first time, to carry out sensitive searches 
for gamma-ray pulsations from all these sources. Even 
before the commissioning phase was complete, Fermi 
already yielded important scientific discoveries in the 
field of pulsar astrophysics, including the detection of 

-http://www.nasa.g0v/mission_pages/GLAST/news/gammaray_best.l 



gamma-ray pulsations from the "Dragonfly" pulsar (PSR 
J202 1-1-3651) 1 22 1, the detection of the first gamma-ray 
ghtch from PSR B 1706^4 [191, and the discovery of 
a new radio-quiet gamma-ray pulsar in the SNR CTA 
1 f3l. With the discovery of a large number of additional 
gamma-ray pulsars 121], Fermi is fulfilling its pre-launch 
expectations by identifying a large number of formerly 
unidentified gamma-ray sources and at the same time 
uncovering a whole new population of pulsars never 
before accessible by any other means. 
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